X- entanglement of PDC photon pairs 
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We investigate the spatio-temporal structure of the bi-photon entanglement in parametric down- 
conversion (PDC) and we demonstrate its non-factorable X-shaped geometry. Such a structure gives 
access to the ultra-broad bandwidth of PDC, and can be exploited to achieve a bi-photon temporal 
localization in the femtosecond range. This extreme localization is connected to our ability to resolve 
the photon positions in the source near-field. The non factorability opens the possibility of tailoring 
the temporal entanglement by acting on the spatial degrees of freedom of twin photons. 
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Parametric down-conversion (PDC) is probably the 
most efficient and widely used source of entangled photon 
pairs, which has been employed in several successful im- 
plementations of quantum communication and informa- 
tion schemes. At the very heart of such technologies lies 
the quantum interference between photonic wave func- 
tions, which depends crucially on the spatio-temporal 
mode structure of the photons. In this work, the issue of 
controlling and tailoring the bi-photon spatio-temporal 
structure is addressed from a peculiar and novel point 
of view, that is, the non factorability in space and time 
of the PDC bi-photon entanglement. The idea comes 
from the context of nonlinear optics, where recent stud- 
ies [l] outlined how in nonlinear media the angular dis- 
persion relations impose a hyperbolic geometry involving 
both temporal and spatial degrees of freedom in a non- 
factorable way. The wave object that captures such a 
geometry is the so-called X-wave (the X being formed by 
the asymptotes of the hyperbola), which is a localized 
and propagation-invariant wave-packet, non separable in 
space and time. The statistical counterpart of the X-wave 
was recently showed to emerge in the X-shaped structure 
of the classical first order coherence function ^2]. 

In this work, we turn our attention to the genuine 
quantum properties of PDC, and we adopt the X-wave 
picture for investigating the spatio-temporal structure of 
the two-photon entangled state. With few exceptions 
[3l H], investigations on the quantum state of PDC have 
been performed to date mostly either in a purely tempo- 
ral [HllllT] or spatial [8l[9l[l0] framework. Our approach, 
based on the non-factorability in space and time of the 
state, will point out relevant elements of novelty, as the 
possibility of tailoring the temporal bandwith of the bi- 
photons by manipulating their spatial degrees of freedom. 
In particular, by resolving their near- field positions, we 
will show that the X-structure opens the access to an 
ultra-broad bandwidth entangled photonic source, with 
a temporal localization in the femtosecond range. Such 
an extreme localization can be applied to increase the 
sensitivity of high precision measurements in the time do- 
main (e.g. in the protocol of clock syncronization [TT] or 
of quantum optical coherence tomography [12 ). Our re- 



sults compare with recent findings reported in [13 , where 
a ~ 7fs Hong-Ou-Mandel dip was observed through the 
use of a quasi-phase-matched nonlinear grating. 

We shall focus on type I PDC, in the low-gain regime 
where single pairs of photons can be detected by coinci- 
dence counts. We remark that the X-structure of entan- 
glement is a general feature of PDC, present also in type 
II and in the high-gain regime [M]. 

The model is basically the same as in [4l, fTO]. A 
quasi-monochromatic and coherent pump field propa- 
gates along the direction z inside a slab of nonlinear x^^^ 
crystal of length /c- ^, ^s{^i z) denote the en- 

velope operators of the pump and the down converted 
signal field, of central frequencies Up and ujs = 
respectively. Here x = {x^y) labels the transverse coor- 
dinates, while t is time. We next pass to the Fourier 
domain: Ai{q,u,z) = J ^ J ^Ai{x,t, z)e-'^-^+'^\ 
i = s^p and extract the fast variation due to their linear 
propagation along the crystal slab: 



Ai{q,u;,z) ■ 
where kiz{q,u) = ^/kf{q,^J) 
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q'^ is the z-component of 
the wave vector of the i-th field, ki{q^u;) being the wave 
number at frequency cj, which for the extraordinary wave 
depends also on the propagation direction, identified by 
the transverse wave- vector q. The fields defined in this 
way have a slow variation along the crystal, arising only 
from the nonlinear interaction. In the low-gain regime we 
can assume that the pump is undepleted by the nonlin- 
ear interaction, so that j^dp{q,u;, z) = 0, and substitute 
its field operator by a c- number field ap{q^u;). In this 
way, the pump evolution along the crystal is described by 
Ap{qjCO^z) = e^^p^^^'^^^ap{q^ij). For the down-converted 
signal field, its propagation along the crystal is described 
by the equation [4 : 
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where g is the dimensionless parametric gain, propor- 
tional to the second-order x^^^ susceptibility, to the crys- 
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tal length Ic and to the pump peak value (the pump 
field has been normalized to its peak value). The phase- 
mismatch function 

(3) 

determines how efficiently a pump photon with trans- 
verse wave- vector q ^ q' and frequency uOp ^ uo ^ uo' is 
down-converted into a pair of photons with transverse 
wave- vectors q^q' and frequencies uOs + Cc;, uOs uo': the 
smaller its modulus, the higher the probability that such 
an elementary process occurs. 

In most experiments in the low gain regime, the quan- 
tity of primary interest is the two-photon correlation, also 
called bi-photon amplitude. We shall study this quantity 
in the spatio-temporal domain, in a plane at the output 
face of the crystal {near-field region)^ that is, we focus on 



ilj{x,t,x',t') = {As{x,t,lc)As{x',t',lc)) . 
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In the low gain limit considered in this work its square 
modulus x', is proportional to the two pho- 

ton coincidence rate G^'^^ (x, t, x', t), which gives the joint 
probability distribution of finding two photons in posi- 
tion X at time t and position x' at time t' ^ respectively. 

For small gains {g <C 1), the propagation equation ([2| 
can be solved perturbatively up to first order in ^, obtain- 
ing the following expression for the bi-photon amplitude 
in the Fourier domain: 

{As{qi,Ui,lc)As{q2,(^2,lc)) = gAp{qi +^2,^i +^2,^c) 



X (2^)-3/2e^ 
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In the literature the same quantity is usually derived 
through a perturbative evaluation of the two-photon 
state vector (see e.g. [3]). 

In order to simplify our results, we consider the limit of a 
nearly plane- wave and monochromatic pump, i.e, a pump 
of duration Tp and waist Wp large enough, so that the de- 
pendence of the phase-mismatch A on gi +g2 and ui -{-002 
(the pump variables) can be neglected. It can be shown 
[I4j that such an approximation holds when Wp and Tp 
are much larger than the spatial walk-off and the tem- 
poral delay due to group velocity mismatch, respectively, 
experienced by the signal and the pump after crossing 
the crystal. Typical values are ~ 300 /xm and ~ 2ps, as 
in the example of a 4 mm /3-barium-borate (BBO) crys- 
tal cut for degenerate type I PDC at 352 nm. Provided 
that we are in such conditions, the bi-photon amplitude 
Q at the crystal output face takes the factorized form: 



^l){x,t,x' ,t') =Ap ( ) tl^py,{x - x\t -t') 
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Apw(<f,cc;) = ksz{q,(^) ^ ksz{-q,-uj) - kp{0,0) {7c) 

is the PWP result for the field correlation function. The 
pump beam profile Ap{x^t^lc) acts thus as a slow mod- 
ulation over the PWP correlation ?/;pw(x, t), as it can be 
expected in the nearly stationary and homogeneous con- 
ditions considered here. 

A first qualitative insight into the problem can be ob- 
tained by considering the usual quadratic expansion of 
the phase matching function, equivalent to adopting the 
paraxial and quadratic dispersion approximations. In the 
case of e-00 phase matching, it takes the form [4j: 



Ap^{q,u)lc - Ao 



(8) 



where Aq = {2ks — kp)lc is the collinear phase mismatch 
at degeneracy, VLq = ^/lJk^, qo = y/ksjlc^ and we used 
the short-hand notation ks = /cs(0,0), kp = /cp(0,0), 
kg = d'^kg/dco'^lo^o. If we extend the validity of such 
an approximation to the entire {q^uj) domain, and we 
use the identity e*^/^sinc(p/2) = dse^^^ ^ the bi-photon 
amplitude ?/^pw(^?, t) can be recasted in the integral form: 



V^pw(r,^) =9 
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where r indicates the radial coordinate. This expression 
clearly evidences the hyperbolic geometry of ?/^pw(r, t): 
the function is indeed constant on the rotational hyper- 
boloids where the argument 
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assumes constant values. However, it can be easily shown 
that ip^^{H) goes as 1/ ^J\H\ for \H\ 0, that is, when 
approaching the asymptotes of the X-structure, where 
(r, t) = q^r'^ — Vt^t^ = 0. This singularity arises from 
the unphysical assumption that the approximation (|8| is 
valid everywhere. 

In order to obtain quantitative results we need there- 
fore to drop the approximation ([8|, and to go beyond 
the paraxial and quadratic dispersion approximations in 
the evaluation of A(g, cj). '0pw(^, t), defined by Eq.([7|), is 
hence numerically calculated by using the complete Sell- 
meier relations [15^ for the refractive indexes. An exam- 
ple of our results is shown by Fig {T] for the case of a type 
I BBO crystal. Since the signal is an ordinary wave, the 
spatial radial symmetry of the problem can be exploited 
to compute Eq.Q by means of a Fourier-Hankel trans- 
form. Correspondingly, the bi-photon amplitude tjjp^ de- 
(6) pends only on the radial coordinate r = \x\. A cut of this 
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FIG. 1: (Color online) (a) and (b): Modulus of the bi-photon 
amplitude in the plane (x,t), clearly displaying its X-shaped 
geometry (the whole 3-D plot has a rotational symmetry in 
space), (c) and (d): Cuts of the coincidence rate |'0pw(^, t)|^ 
along the temporal (c) and the transverse coordinate axis (d) 
[indicated by the two dashed lines in frame (b)]. The width of 
the peaks shows the relative temporal and spatial localization 
of bi-photons. BBO crystal, cut at 33.436° for type I PDC; 
g — 10~^, \p — 352nm, and Ic — 4mm. 



quantity in the (x, t) plane is displayed by frames (a) and 
(b) of Figjl] (the whole 3-dimensional plot have a radial 
symmetry in the space domain, and has therefore a bi- 
conical geometry). A clear X-shaped structure emerges: 
this non-trivial shape of the spatio-temporal two-photon 
correlation, that we shall call X- entanglement cdin be con- 
sidered the counterpart, at the level of quantum noise, of 
the nonlinear X- waves [1 . Similar results hold for a va- 
riety of phase-matching conditions, corresponding to dif- 
ferent crystal orientations, and for a wide range of crystal 
lengths [14 . 

A remarkable characteristic of the X-entanglement is 
the unusually small width of the spatio-temporal corre- 
lation peak, which corresponds to a strong relative local- 
ization of twin photons both in time and space. The two 
lower frames of Fig{l] plot cuts of the two-photon coin- 
cidence rate |?/^pw(^, t)P along the temporal and spatial 
axis, respectively. The spatial localization is remarkable 
but not impressive, as displayed by the spatial profile 
|V^pw(x, 0) p in Fig{ljd), which has a FWHM of - 2.9 /im. 
More impressive, and in a sense, unexpected, is the tem- 
poral relative localization of twin photons, which can 
be appreciated from the temporal profile |?/^pw(0,t)p in 
Figjljc), which is as narrow as 4.4 fs. Such an ultra-short 
two-photon localization emerges spontaneously from a 
nearly monocromatic pump, as a consequence of the 
ultra-broad bandwidth of PDC phase-matching, which in 
principle extends over the optical frequency cOp ^ 5 10^^ 
Hz. Notice that, in order to account for e.g. the finite 



bandwidth of detection, in our calculations we include 
a super-gaussian frequency filter centered at degeneracy. 
The 4.4 fs width of the temporal peak is in practice de- 
termined by the width of this frequency filter, (see Figj2] 
for a comparative view of the frequency bandwidths in- 
volved) . 

It is interesting to compare our results with the typ- 
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FIG. 2: (Color online) Plot of |y(^cc,cj)|, showing the phase 
matching curve. The dashed vertical lines show the frequency 
bandwidth selected by the filter and the arrows indicate the 
different bandwidths involved (same parameters as in Figjl]). 

ical ~ 100 fs temporal localization of the coincidence 
rate measured in the far- field zone, by collecting twin 
photons that propagate at symmetric directions q and 
— within a small angular bandwidth. In that case, 
the measured quantity is proportional to \V{q^t)\'^ = 
I J |^g-*a;^y(^^^^j|2 j|.g temporal width is determined 
by the inverse of the bandwidth of the spectrum V{q^ co) 
at fixed i.e. the narrow (10^^ — lO^^Hz) thickness of 
the curve in Figjg] This bandwidth can be roughly eval- 
uated as rto for q/qo « 1, and as Vt^q^/q for q/q^ > 1. 

— 1/2 

Clearly, since VLq scales as /c , the shorter the crys- 
tal, the stronger the temporal localization in the far- 
field; however, a far-field localization in the femtosecond 
range would require a crystal as short as ~ 50 /im, with a 
strongly reduced down- conversion efficiency. Conversely, 
in our case, the detection of coincidences in the near-field 
gives in principle access to the full (~ lO^^Hz) bandwidth 
of phase matching even for a long crystal. 
It is however important to stress that such an extreme 
temporal localization of twin photons relies on the abil- 
ity to resolve their positions in the near-field plane. In- 
deed, a measurement collecting all the photons over the 
beam cross-section, without discriminating their posi- 
tions, is characterized by the integrated coincidence rate 
J |?/^pw(^, t)P, reproduced by the dashed curve in 
FigjSj which has a width of ~ lOOfs. This may appear 
surprising, because in this measurement all the photons 
at the different frequencies within the phase-matching 
are collected. However, the identity / d?x\ilj^^{x^t)\^ = 
J dPq\V{q^t)\^ ^ shows that in this case the coincidence 
rate takes the form of an incoherent superposition of the 
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FIG. 3: Solid line: |'0pw(^ = 0, t) coincidence rate when 
photon positions are resolved in the near field. Dashed line: 
j dPx |'0pw(^, coincidence rate measured without resolv- 
ing photon positions 
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FIG. 4: (Color online) Effect of spatial filtering on the X- 
entanglement: (a) Temporal correlation peak |'0pw(O,t)|^ in 
the presence of a spatial filter, that cuts the angular spectrum 
at an angle amax- The two insets show the full X-correlation 
for amax — 2*^ (b), and in the absence of any filter (c). 



probabilities of detecting a pair of photons at a given q 
and has therefore the same ~ lOOfs temporal localization 
as the far-field coincidence rate at fixed q. Conversely, 
by resolving the near-field positions of twin photons, 
the measured quantity is |V^pw(0,t)p = | / -^4^V {q,t)\^ , 
which corresponds to a coherent superposition of the 
probability amplitudes at a given q (i.e. at a given fre- 
quency due to the the angle- frequency relation imposed 
by phase matching) , and therefore allows a stronger tem- 
poral localization. 

The non-factorability in space and time of the X- 
entanglement thus opens the relevant possibility of tailor- 
ing the temporal bandwidth of the bi-photons by acting 
on their spatial degrees of freedom. As a more specific ex- 
ample, let us consider the effects of spatial filtering on the 
temporal correlation. Let us assume that a 4/ lens sys- 
tem is employed to image the near-field of the PDC flu- 
orescence, and that a circular aperture of radius is lo- 
cated in the far-field 2/ plane, acting as a filter that cuts 
all the angular spectrum at a > amax = arcsin(ra//). 
Fig. [4] shows the effect of such a spatial filter on the 
temporal correlation peak. While in the absence of any 
spatial filter the correlation shows a strong temporal lo- 



calization, as the angular bandwidth is reduced by spatial 
filtering, the two-photon correlation broadens in time. 
This is a clear effect of the non-factorability of the corre- 
lation, because, thanks to the shape of the angular spec- 
trum shown in Figj2| a spatial filter that cuts the angular 
bandwidth has also the effect of cutting the frequency 
bandwidth. 

In conclusions, this work demonstrates the hyperbolic 
geometry underlying the two-photon PDC entanglement 
and its non-factorability with respect to space and time. 
As for the X- waves encountered in nonlinear optics, the 
X-shape of the bi-photon correlation is imposed by the 
phase-matching mechanism governing the PDC process, 
and following this analogy we coined the name of X- 
entanglement. The key element of novelty brought by 
this structure is its extreme localization, with correla- 
tion times and correlation lengths in the femtosecond and 
micrometer range, respectively. The strong temporal lo- 
calization is determined by the full extent of the PDC 
bandwidth, rather than by the bandwidth ~ fto charac- 
terizing the PDC far-field. For this reason, a near field 
measurement scheme able to resolve spatially the coinci- 
dences would provide a powerful tool for high-precision 
measurements, capable of improving substantially the 
resolution power in the temporal domain with respect 
to standard schemes. Furthermore, we have shown how 
the non-factorability of the structure gives the possibility 
of tailoring the temporal entanglement through manipu- 
lations of the spatial degrees of freedom. 
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